home *** CD-ROM | disk | FTP | other *** search
/ Fritz: All Fritz / All Fritz.zip / All Fritz / FILES / PROGMISC / PCSSP.LZH / PC-SSP.ZIP / STATMISC.ZIP / BISER.FOR next >
Text File  |  1985-12-31  |  3KB  |  131 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C     SUBROUTINE BISER
  5. C
  6. C     PURPOSE
  7. C        TO COMPUTE THE BISERIAL CORRELATION COEFFICIENT BETWEEN TWO
  8. C        CONTINUOUS VARIABLES WHEN ONE OF THEM HAS BEEN ARTIFICIALLY
  9. C        DICHOTOMIZED.
  10. C
  11. C     USAGE
  12. C        CALL BISER (N,A,B,HI,ANS,IER)
  13. C
  14. C     DESCRIPTION OF PARAMETERS
  15. C        N    - NUMBER OF OBSERVATIONS
  16. C        A    - INPUT VECTOR OF LENGTH N CONTAINING THE CONTINUOUS
  17. C          VARIABLE
  18. C        B    - INPUT VECTOR OF LENGTH N CONTAINING THE DICHOTOMIZED
  19. C          VARIABLE
  20. C        HI    - INPUT - NUMERICAL CODE TO INDICATE THE HIGHER CATEGORY
  21. C          OF THE DICHOTOMIZED VARIABLE.  ANY VALUE IN VECTOR B
  22. C          EQUAL TO OR GREATER THAN HI WILL BE CLASSIFIED INTO
  23. C          THE HIGHER CATEGORY.
  24. C        ANS - OUTPUT VECTOR OF LENGTH 8 CONTAINING THE FOLLOWING
  25. C          ANS(1) - MEAN OF VARIABLE A
  26. C          ANS(2) - STANDARD DEVIATION OF VARIABLE A
  27. C          ANS(3) - PROPORTION OF THE CASES IN THE HIGHER
  28. C               CATEGORY OF VARIABLE B
  29. C          ANS(4) - PROPORTION OF THE CASES IN THE LOWER
  30. C               CATEGORY OF VARIABLE B
  31. C          ANS(5) - MEAN OF VARIABLE A FOR THOSE CASES FALLING
  32. C               INTO THE HIGHER CATEGORY OF VARIABLE B
  33. C          ANS(6) - MEAN OF VARIABLE A FOR THOSE CASES FALLING
  34. C               INTO THE LOWER CATEGORY OF VARIABLE B
  35. C          ANS(7) - BISERIAL CORRELATION COEFFICIENT
  36. C          ANS(8) - STANDARD ERROR OF BISERIAL CORRELATION
  37. C               COEFFICIENT
  38. C        IER -  1, IF NO CASES ARE IN THE LOWER CATEGORY OF VARIABLE
  39. C          B.
  40. C          -1, IF ALL CASES ARE IN THE LOWER CATEGORY OF
  41. C          VARIABLE B.
  42. C          0, OTHERWISE.
  43. C          IF IER IS NON-ZERO, ANS(I)=10**75,I=5,...,8.
  44. C
  45. C     REMARKS
  46. C        THE VALUES OF THE DICHOTOMIZED VARIABLE, B, MUST BE IN
  47. C        NUMERIC FORM.  THEY CANNOR BE SPECIFIED BY MEANS OF
  48. C        ALPHABETIC OR SPECIAL CHARACTERS.
  49. C
  50. C     SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  51. C        NDTRI
  52. C
  53. C     METHOD
  54. C        REFER TO P. HORST, 'PSYCHOLOGICAL MEASUREMENT AND
  55. C        PREDICTION', P.95-96 (WADSWORTH, 1966).
  56. C
  57. C     ..................................................................
  58. C
  59.       SUBROUTINE BISER (N,A,B,HI,ANS,IER)
  60. C
  61.       DIMENSION A(1),B(1),ANS(1)
  62. C
  63. C     COMPUTE MEAN AND STANDARD DEVIATION OF VARIABLE A
  64. C
  65.       IER=0
  66.       SUM=0.0
  67.       SUM2=0.0
  68.       DO 10 I=1,N
  69.       SUM=SUM+A(I)
  70.    10 SUM2=SUM2+A(I)*A(I)
  71.       FN=N
  72.       ANS(1)=SUM/FN
  73.       ANS(2)=(SUM2-ANS(1)*SUM)/(FN-1.0)
  74.       ANS(2)= SQRT(ANS(2))
  75. C
  76. C     FIND PROPORTIONS OF CASES IN THE HIGHER AND LOWER CATEGORIES
  77. C
  78.       P=0.0
  79.       SUM=0.0
  80.       SUM2=0.0
  81.       DO 30 I=1,N
  82.       IF(B(I)-HI) 20, 25, 25
  83.    20 SUM2=SUM2+A(I)
  84.       GO TO 30
  85.    25 P=P+1.0
  86.       SUM=SUM+A(I)
  87.    30 CONTINUE
  88.       ANS(4)=1.0
  89.       ANS(3)=0.0
  90.       Q=FN-P
  91.       IF (P) 35,35,40
  92.    35 IER=-1
  93.       GO TO 50
  94.    40 ANS(5)=SUM/P
  95.       IF (Q) 45,45,60
  96.    45 IER=1
  97.       ANS(4)=0.0
  98.       ANS(3)=1.0
  99.    50 DO 55 I=5,8
  100.    55 ANS(I)=1.E38
  101.       GO TO 65
  102.    60 ANS(6)=SUM2/Q
  103.       P=P/FN
  104.       Q=1.0-P
  105. C
  106. C     FIND ORDINATE OF THE NORMAL DISTRIBUTION CURVE AT THE POINT OF
  107. C     DIVISION BETWEEN SEGMENTS CONTAINING P AND Q PROPORTIONS
  108. C
  109.       CALL NDTRI (Q,X,Y,ER)
  110. C
  111. C     COMPUTE THE BISERIAL COEFFICIENT OF CORRELATION
  112. C
  113.       R=((ANS(5)-ANS(1))/ANS(2))*(P/Y)
  114. C
  115. C     COMPUTE THE STANDARD ERROR OF R
  116. C
  117.       ANS(8)=( SQRT(P*Q)/Y-R*R)/SQRT(FN)
  118. C
  119. C     STORE RESULTS
  120. C
  121.       ANS(3)=P
  122.       ANS(4)=Q
  123.       ANS(7)=R
  124. C
  125.    65 RETURN
  126.       END
  127.  ANS(6)=SUM2/Q
  128.       P=P/FN
  129.       Q=1.0-P
  130. C
  131.